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Abstract 

Anticausal  inversion  of  HR  transfer  functions  has 
gained  importance  in  recent  years,  in  the  efficient  im¬ 
plementation  of  HR  digital  filter  banks.  In  this  paper 
we  first  introduce  the  idea  of  a  causal  dual,  as  an  inter¬ 
mediate  step  in  the  implementation  of  anticausal  IIR 
inverses.  With  time  reversal  operators  at  the  input 
and  output  of  the  casusal  dual,  we  get  the  anticausal 
inverse  of  the  original  structure.  The  causal  dual  elim¬ 
inates  the  need  for  similarity  transformations,  during  a 
key  step  called  blockwise  state  transfer,  in  implement¬ 
ing  anticausal  inverses.  In  the  paper  we  identify  effi¬ 
cient  structures  for  causal  duals  of  standard  structures 
like  the  direct-form,  cascade-form,  coupled  form,  and 
IIR  lattice  structures,  including  the  tapped  lattice. 

1.  INTRODUCTION 

Anticausal  inversion  of  IIR  transfer  functions  has  gained 
importance  in  recent  years,  in  the  implementation  of 
digital  filter  banks  [l]-[6].  This  is  due  to  a  certain 
class  of  IIR  digital  filter  banks,  in  which  the  synthesis 
filters  are  required  to  be  anticausal  if  a  stable  recon¬ 
struction  system  (synthesis  bank)  is  desired.  Thus, 
Fig.  1.1  shows  the  polyphase  form  [7]  of  a  commonly 
used  two  channel  filter  bank.  This  nas  perfect  recon¬ 
struction  [x(n)  =  x(n)]  if  and  only  if  J?i(z)  =  lfEi{z). 
If  Ei{z)  does  not  have  all  zeros  inside  the  unit  circle, 
then  the  causal  impulse  response  of  Ri{z)  is  unstable. 
In  some  cases  (e.g.,  when  Ei{z)  is  allpass  with  poles 
inside  the  unit  circle),  1/Ei{z)  has  all  poles  outside, 
so  that  an  anticausal  impulse  response  of  1/Ei{z)  is 
stable.  This  explains  the  interest  in  implementing  an¬ 
ticausal  inverses.  This  has  been  shown  to  be  possible 
[3]-[6],  even  in  real  time  with  infinite  duration  inputs, 
as  long  as  we  perform  the  processing  in  blocks  of  length 
L,  and  properly  transmit  the  state  variables  from  Ei(z) 
to  1/Ei{z)  at  the  end  of  each  block. 

Consider  a  causal  stable  ATth  order  filter  G{z),  im¬ 
plemented  using  a  minimal  structure  (direct-form,  lat¬ 
tice,  etc.),  with  state  space  description 
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where  x(n)  =  [xi(n)  ...  xjv(n)  is  the  state  vec¬ 
tor  [xi(n)  are  the  outputs  of  delays],  it(n)  is  the  filter 
input  and  y{n)  the  filter  output.  The  transfer  function 
is  given  by  G(z)  =  D  +  C(zl  —  A)“^B  and  has  an 
anticausal  inverse  [i.e.,  there  is  an  anticausal  inverse 
z-transform  of  1/G(z)]  if  and  only  if  the  realization 
maim  72.  is  nonsingular  [5].  We  assume  this.  Consider 
the  following  “dual”  casual  state  space  equations 
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where  72.  =  72“^.  Imagine  we  run  the  original  sys¬ 
tem  (1.1)  for  a  duration  of  L  samples,  with  input 
u(0) . . .  u{L  —  1)  and  get  the  output  j/(0) . . .  y{L  -  1) 
and  final  state  x(L).  If  we  run  the  state  space  equations 
(1.2)  by  setting  the  initial  state  to  be  x(£)  =  x(I-), 
and  the  input  sequence  to  be  the  time  reversed  block 
y{L  —  1)  •  •  •3/(0),  we  can  show  that  the  output  of  this 
system  for  the  duration  of  this  block  is  u{L—l) . . .  u(0). 
In  this  way  we  can  invert  (1.1),  i.e.,  invert  G(z).  Since 
the  time  reversal  is  blockwise,  it  can  be  realized  in 
practice  with  finite  latency.  By  repeating  the  above 
process  block  by  block,  we  can  implement  the  anti¬ 
causal  inverse  1/G{z)  of  the  original  filter  G{z),  even 
for  infinitely  long  inputs. 

It  can  be  shown  that  the  transfer  fimction  of  the 
causal  system  (1.2)  is  1/G{z~^).  We  can  schematically 
express  the  desired  anticausal  inverse  1/G{z)  as  in  Fig. 
1.2,  in  terms  of  ideal  time  reversal  (TR)  operators.  In 
a  practical  implementation  of  the  “anticausal”  inverse 
1/G{z),  the  time  reversals  are  done  blockwise,  and  the 
causal  system  1/G{z~^)  implemented  using  (1.2).  The 
initial  condition  for  the  state  recursion  (1.2)  at  the  be¬ 
ginning  of  each  block  is  taken  as  the  final  state  from 
the  previous  block  of  the  state  recursion  (1.1)  for  G{z). 
This  state  transfer  is  crucial. 

Aim  of  the  Paper 

Note  that  \/G[z~^)  can  be  realized  using  any  structure 
(direct-form,  lattice,  etc.).  The  state  space  description 

of  these  structures  are  related  to  (A,  B,  C,  D)  by  sim¬ 
ilarity  transforms.  If  an  arbitrary  structure  is  used 
for  1/G(z~^),  then  the  state  transfer  described  above 
should  be  done  only  after  a  similarity  transformation 
at  the  beginning  of  each  block. 


Causal  dual  of  a  structure.  A  structure  for 
1/G(z“^)  which  has  precisely  the  state  space  descrip¬ 
tion  (A,  B,  C,  D)  will  be  called  a  causal  dual  of  the 
structure  (A,B,  C,jD)  for  G{z).  If  the  causal  dual  is 
used  for  implementing  1/G{z~^),  no  similarity  trans¬ 
formations  are  necessary  during  state  transfer.  In  this 
paper  we  will  identify  efficient  structures  for  the  causal 
duals  of  standard  structures,  such  as  the  direct-form, 
cascade- form,  coupled  form,  and  the  family  of  lattice 
structures,  including  the  tapped  lattice. 

2.  CAUSAL  DUAL  FOR  THE  DIRECT-FORM 

Fig.  2.1(a)  shows  the  direct  form  for  an  iVth  order  filter 
G(i).  Prom  [5]  we  know  that  the  inverse  filter  1/G(z) 
has  an  anticausal  impulse  response  (more  simply,  G(z) 
has  an  anticausal  inverse)  if  and  only  if  ^  0;  we 
therefore  assume  this.  With  the  state  variables  Xi{n) 
as  indicated  in  the  figure,  the  matrix  A  is  in  compan¬ 
ion  form  (e.g.,  see  [7]).  Writing  down  the  realization 
matrix  TZ,  we  can  explicitly  invert  it,  and  verify  that 
the  structure  in  Fig.  2.1(b)  has  the  realization  matrix 
'Rr^,  and  is  therefore  the  causal  dual.  Thus  the  causal 
dual  can  be  obtained  by  making  a  simple  set  of  changes 
to  the  multipliers,  and  renumbering  the  state  variables 
in  reverse  order.  For  the  special  case  of  allpass  filters 
the  causal  dual  is  even  simpler,  since  p„  =  For 

real  coefficient  allpass  filters,  the  direct  form  structure 
is  its  own  causal  causal  dual  —  only  the  states  need  to 
be  renumbered.  This  is  consistent  with  the  anticausal 
inversion  of  direct-form  reported  in  [8]. 

Cascade-form  structures.  The  cascade-form, 
which  is  a  cascade  connection  of  second  order  direct 
form  structures,  is  well-known  for  its  simplicity  and 
relatively  better  robustness  to  quantization  [10].  More 
generally,  let  Gi{z)  and  Gaiz)  be  causal  transfer  func¬ 
tions  in  cascade  (Fig.  2.2),  implemented  with  struc¬ 
tures  having  realization  matrices  Tlx  and  7^2-  Then 
from  the  block  of  L  outputs  y2{n) . .  .y2{n  +  L  —  1) 
and  the  state  X2(n  4-  L)  of  G2(z),  we  can  recover  the 
block  of  L  inputs  yi{n) . . . yx{n  +  L  —  1)  of  the  system 
G2(z).  We  can  then  use  yi{n) ... yx{n  +  L  —  1)  (out¬ 
put  block  of  Gi{z))  and  the  state  vector  xi(n  4-  L)  of 
the  system  Gi(z)  to  recover  the  primary  input  block 
u(n) . . .  u(n  +  L  —  1).  This  is  equivalent  to  connecting 
the  anticausal  inverses  in  reverse  order. 

In  fact,  it  can  be  shown  that  the  causal  dual  of  a 
cascade  is  simply  the  cascade  of  the  individual  causal 
duals,  in  reverse  order.  If  we  implement  this  cascaded 
causal  dual  structure  with  blocWise  time  reversal  of 
its  primary  input  and  output,  we  obtain  the  anticausal 
implementation  of  the  complete  system  1/G{z).  Thus, 
regardless  of  the  number  of  filters  connected  in  cascade, 
time  reversal  is  necessary  only  at  the  primary  input 
and  output  nodes,  and  not  at  the  intermediate  inputs 
and  outputs. 

3.  CAUSAL  DUAL  FOR  THE  COUPLED-FORM 

Fig.  3.1(a)  shows  the  coupled  form  structure,  whose 
robustness  to  quantization  is  well-known  [10].  The 
poles  of  this  system  are  at  pe^^  and  With  the 


output  node  y(n)  as  indicated,  the  transfer  function 
is  G{z)  =  psm6z~‘^ /{I  —  2pcos9z~^  -f  p'^z~^).  Upon 
analysis  we  find  that  the  realization  matrix  is 
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Fig.  3.1(b)  shows  the  structure  with  this  realization 
matrix  (where  m  =  psin0-f  pcot^cos^).  This  is  the 

causal  dual  of  Fig.  3.1(a).  Note  that  A  (the  top-left 
2x2  matrix  above)  has  all  eigenvalues  =  0.  This  is 
because  1/G(^;“^)  is  FIR. 

The  Tapped  Coupled  Form.  Fig.  3.2(a)  shows 
the  tapped  coupled  form.  This  has  the  extra  multipli¬ 
ers  (tap  coefficients)  D,  ci  and  C2  which  can  be  used 
to  obtain  arbitrary  numerators.  The  realization  matrix 
now  becomes 
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The  inverse  of  this  matrix  is 
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which  can  be  rewritten  as 
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with  a  =  (ci  —  £>pcos0)/psin0,  13  =  C2  +  Dpsin^  — 
ap  cos  6.  Thus  TZf^p  can  be  written  as 
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(3.3) 


The  structure  with  this  realization  matrix  can  be  ob¬ 
tained  from  Fig.  3.1(b)  by  replacing  the  input  tt(n) 
with  —{D/0)xi{n)  —  {a/l3)x2{n)  +  {l//3)u{n),  as  shown 
in  Fig.  3.2(b).  This  is  therefore  the  causal  dual  of  the 
tapped  coupled-form. 


4.  CAUSAL  DUALS  FOR  LATTICE  STRUCTURES 

In  the  case  of  the  direct  form  structure,  we  could  obtain 
the  causal  dual  simply  by  inspection  (Fig.  2.1).  For 
arbitrary  structures,  the  rules  could  be  more  compli¬ 
cated,  as  shown  by  the  coupled  form  example.  We  will 
see  next  that  for  the  lattice  structures,  which  are  well- 
known  for  many  good  robustness  properties  [11]-[12], 
the  causal  dual  can  be  obtained  again  by  inspection! 
The  lattice  structure  is  shown  in  Fig.  4.1(a).  Here 
Gn{z)  =  Y{z)/U{z)  is  a  causal  Nth  order  allpass  filter. 
The  rectangular  building  block  labelled  km  can  take 
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several  possible  forms  (Fig.  4.1(b)).  While  the  nor¬ 
malized  structure  offers  automatic  scaling  and  normal¬ 
ization  of  noise  gains  [11],  the  two- multiplier  version 
offers  economy.  For  a  given  set  of  lattice  coefficietns 
{km},  the  allpass  function  Gn{z)  is  the  same  regardless 
of  which  of  these  building  blocks  is  used.  The  lattice 
coefficients  km  are  possibly  complex,  but  |fcm|  <  1,  so 

Gn{z)  is  stable.  The  multiplier  km  —  \/l  —  [AimP  is 
always  real. 

With  the  state  variables  Xm{n),  the  input  u{n)  and 
output  y{n)  indicated  as  shown  in  Fig.  4.1(a),  let  the 
realization  matrix  of  this  Nth.  order  structure  be  de¬ 
noted  as  TZn-  We  will  show  that  the  causal  dual  (i.e., 
the  structure  with  realization  matrix  TZJ;^)  is  as  shown 
in  Fig.  4.2,  where  the  boxes  k^  are  the  bailing  blocks 
in  the  original  structure  Fig.  4.1,  but  with  coefficients 
conjugated.  Thus  the  causal  dual  structure  is  obtained 
from  the  original  structure  by  conjugating  the  multipli¬ 
ers,  and  moving  the  delays  from  the  bottom  rails  to  the 
top  rails. 

For  Fig.  4.2,  define  the  state  vector  x(n)  and  the 
state  space  description  (A,B,  C,£)),  and  let  the  real¬ 
ization  matrix  be  denoted  TZ^.  First  consider  the  nor¬ 
malized  lattice.  In  this  case  TZn  is  unitary  [13],  so  we 

only  have  to  show  TZn  ■=  7?.^  (transpose  conjugate). 
This  is  based  on: 

Lemma  3.1.  The  (AT  +  1)  x  (iV  -f- 1)  realization 
matrix  IZn  for  the  lattice  of  Fig.  4.1(a)  with  normal¬ 
ized  building  blocks  can  be  expressed  as  a  product  of 
N  unitary  matrices: 


Real  coefficient  case.  If  Gn{z)  has  real  coeffi¬ 
cients,  then  km  are  real.  So  the  causal  dual  is  identical 
to  the  original  structure,  except  that  delays  have  to  be 
moved  up. 

Denormalized  lattice  sections.  Suppose  the  lat¬ 
tice  structure  of  Fig.  4.1(a)  uses  denormalized  sections 
(e.g.,  the  two  multiplier  section  of  Fig.  4.1(b),  or  one- 
multiplier  sections  which  are  available  for  real  coeffi¬ 
cient  case  [7]).  Then  the  structure  is  related  to  the 
normalized  lattice  via  a  diagonal  similarity  transfor¬ 
mation.  Using  this,  it  is  possible  to  show  that  if  Fig. 
4.1(a)  represents  a  denormalized  lattice,  then  Fig.  4.2 
still  represents  its  causal  dual.  That  is,  we  simply  con¬ 
jugate  coefficients  and  move  the  delays  to  upper  rails. 

The  Tapped  Lattice  Structure.  The  tapped  lat¬ 
tice,  which  can  be  used  to  realize  an  arbitrary  HR 
transfer  function  H{z)  is  shown  in  Fig.  4.3(a).  Here 
H {z)  has  denominator  equal  to  the  that  of  the  allpass 
filter  Gn{z).  The  tap  coefficients  can  be  chosen  to 
realize  the  numerator  of  H{z).  We  will  show  that  the 
causal  dual  of  this  is  given  by  Fig.  4.3(b).  (This  as¬ 
sumes  ajv+i  ^  0.  If  ccAr+i  =  0,  the  numerator  of  H{z) 
has  a  smaller  order  than  the  denominator,  and  an  an- 
ticausal  inverse  does  not  exist  [5].  The  causal  dual 
would  then  be  of  no  interest.) 

For  fixed  input  u{n),  the  state  variables  Xi(n)  in 
Figs.  4.3(a)  and  4.1(a)  are  identical.  In  Fig.  4.3(a) 

yi(«)  =  I^ili  Q:iXi(n  -fl)  -t-  Q:;v-f-i2/(n)-  Thus  the  real¬ 
ization  matrix  Ttorb  for  tne  arbitrary  transfer  function 
H{z)  is  related  to  the  realization  matrix  TZaii  of  the 
allpass  filter  G/^(z)  as 
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where  the  unitary  matrix 
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Proof.  It  can  be  shown  that  realization  matrix  of 
the  m-stage  lattice  and  the  (m  —  l)-stage  lattice  are 
related  as 
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For  m  =  1  we  can  explicitly  verify  that  7Zi  =  0i.  The 
lemma  follows  from  this  by  induction.  V  VV 

Theorem  3.1.  Causal  dual  for  normalized  lat¬ 
tice.  The  causal  dual  of  the  normalized  lattice  struc¬ 
ture  (Fig.  4.1(a),  with  normalized  building  blocks)  is 
given  by  the  lattice  of  Fig.  4.2.  ^ 

Proof.  We  can  show  that  the  realization  matrices 
for  the  m-stage  lattice  and  (m  —  l)'Stage  lattice  in  Fig. 
4.2  are  related  as 
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and  we  can  explicitly  verify  that  TZi  =  ©{.  From  these 
we  can  obtain  Hn  =  =  TZj/' .  V  VV 
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where  a  =  [  ai  02  ...  ajv  ] .  Thus 
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Since  TZ^i]  corresponds  to  Fig.  4.2,  the  above  realiza¬ 
tion  matrix  corresponds  to  the  structure  Fig.  4.3(b), 
which  is  therefore  the  causal  dual  of  Fig.  4.3(a).  This 
has  taps  in  a  feedback  loop.  Whether  this  is  stable 
depends  on  the  numerator  of  H{z). 
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Fig.  1.1.  A  two-channel  filter  bank  in  polyphase  form. 
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Fig.  1.2.  The  anticausal  inverse  represented  in  terras 
of  a  causal  system  and  time-reversal  operators. 


ytn) 


Fig.  3.1.  (a)  The  coupled-form  structure,  and  (b)  its  causal  dual. 


Fig.  3.2.  (a)  The  ts^ped  coupled  form,  and  (b)  its  causal  dual. 


normalized  block  2-multiplier  block 

Fig.  4.1.  (a)  The  general  form  of  HR  allpass  lattice,  and 
(b)  two  possible  building  blocks. 


Fig.  4.2.  The  causal  dual  of  the  HR  lattice  structure. 


Fig.  2.1  (a)  The  direct-form  structure,  and  (b)  its  causal  dual. 


Fig.  2.2.  A  cascaded  system. 


Fig.  4.3.  (a)  The  tapped  IIR  lattice  and  (b)  its  causal  dual. 


4 


